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Abstract 



The K2K long-baseline neutrino oscillation experiment uses a Scintillating Fiber 
Detector (SciFi) to reconstruct charged particles produced in neutrino interactions 
in the near detector. We describe the track reconstruction algorithm and the per- 
formance of the SciFi after three years of operation. 
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1 K2K Experiment 



The K2K experiment is the first long-basefine neutrino oscillation experiment 
in operation. Muon neutrinos (i^^) are produced via the decay of secondary 
pions from interactions of the KEK 12-GeV proton beam in an aluminum 
target. The neutrino beam is directed through a near detector system at KEK 
to the Super-Kamiokande underground detector[l] 250km away. By comparing 
the neutrino fluxes and energy spectra measured at the near detector and at 
Super-Kamiokande, we are able to investigate neutrino oscillations. Details of 
the K2K experiment have been presented elsewhere [2]. 

The near detector, located 300 m downstream of the pion production point, 
consists of a Ikt water Cherenkov detector, a scintillating flber tracking detec- 
tor (SciFi), a hodoscope of plastic scintillation counters (VETO) around the 
SciFi, a lead glass detector (LG), and a muon range detector (MRD) [3]. It 
is illustrated in Fig. 1. The main purpose of the near detector is to measure 
the profile, fiux, and energy spectrum of the neutrino beam. Among the near 
detector components, the SciFi can provide tracking capability good enough 
to discriminate neutrino interaction types and to measure energy and cross 
sections of the muon neutrino. 

The K2K experiment started taking data in April, 1999. This paper will de- 
scribe the performance of the Scintillating Fiber Detector after three years 
of operation. A description of the K2K scintillating fiber detector is given in 
Section 2. Hit rcconstrTiction algorithm and performance are described in Sec- 
tion 3. Track and vertex reconstruction is described in Section 4. Performance 
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Fig. 1. The layout of the K2K near detector. It consists of a 1-kt water Cherenkov 
detector, a scintillating fiber detector, a hodoscope of plastic scintillation counters, 
a lead glass detector, and a muon range detector. 

of the neutrino event reconstruction algorithms are presented in Section 5. We 
conclude in Section 6. 



2 Scintillating Fiber Detector 

The SciFi is a charged particle tracking detector using scintillating fiber tech- 
nology. Scintillating fiber detectors have been used for charged particle track- 
ers by other experiments [4]. A scintillating fiber detector was chosen for this 
application since it was necessary to have high tracking efficiency for recon- 
structing charged particles and their vertex positions when more than two 
visible particles are produced in neutrino interactions. Using the detector, we 
can separate quasi-elastic events from inelastic events, and define the fiducial 
volume of the neutrino interactions with 1% accuracy. 

The water is contained in a thin aluminium structure, and the water layers 
are interleaved with scintillating fiber layers. Water is chosen as a target ma- 
terial in order to minimize systematic uncertainties in neutrino cross-section 
when the event rates are compared between the SciFi and Super-Kamiokande 
detectors. A scintillating fiber layer provides the hit-position measurement in 
both horizontal and vertical directions [^J An individual scintillation sheet has 

^ In the K2K coordinate system, the z-axis is along the nominal neutrino beam di- 
rection, y-axis is the vertical direction pointing upwards, and x-axis is the horizontal 
direction. 
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Fig. 2. A schematic view of the K2K scintillating fiber detector. The detector has 
the dimension of 260cmx260cmx 173cm in (x, y, z). The central fiducial mass region 
is defined to be 220cmx 220cm x 159cm in {x,y,z). 

a thickness of 1.3mm, consisting of two staggered layers of scintillating fibers. 
Each fiber has a diameter of 0.7mm. The vertical and horizonal fiber sheets 
are separated by a 1.6-cm thick honeycomb panel. There are 20 scintillating 
fiber layers in total, and the adjacent layers are spaced 9cm apart along the 
neutrino beam direction. 

One layer of 240cmx240cmx6cm water target, consisting of 15 aluminium 
tanks, is interleaved with two fiber layers. The wall thickness of the aluminium 
tank is 1.8mm. The weights of water, aluminium, and the honeycomb panel 
are 6.0tons, 1.4tons, and O.Stons, respectively. A schematic view of the SciFi 
is shown in Fig. 2. 

For processing light signals, a group of fibers are bundled and coupled to an 
image intensifier tube (IIT), an optical lens, and a CCD camera as shown 
in Fig. 2. A total of 24 IIT's are used to read out all 274080 fibers in the 
SciFi. In order to find the correspondence between the hit fiber and the CCD 
images, a certain number of selected fibers are illuminated periodically at the 
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Fig. 3. (a) A typical CCD image of an IIT from a cosmic-ray event. Small distortions 
of the CCD image are introduced due to the position dependence of pixel size at 
the IIT-CCD contact, (b) An expanded view of (a) for a typical SciFi hit. The size 
of hit pixel is proportional to brightness in the ADC. Hit fibers are represented by 
the thick circles. A "SciFi hit" is defined as a cluster of neighboring hit fibers. 

opposite side of IIT's by an electro-luminescent (EL) plate. This is referred 
to as the EL calibration. A more detailed description of the SciFi detector is 
found elsewhere [5]. 



3 Hit reconstruction 

3.1 Hit Finding Algorithm 

The CCD pixel images need to be handled to reconstruct hits in the SciFi. 
Fig. 3(a) shows a typical CCD image of an IIT from a cosmic-ray event. A 
fiber bundle, consisting of ten fiber sheets, is arranged to cover a sensitive 
region of an IIT effectively. Fig. 3(b) is an expanded view of (a), illustrating 
a typical "SciFi hit", when a cosmic-ray muon passes through the fiber sheet. 
The hit reconstruction is performed in the following four steps: 

1. Loading Hit Pixels The raw data of the SciFi consist of hit pixels in 
the CCD coordinates {x, y) and 8 bits of ADC brightness. These are all 
encoded in one byte per pixel. Typically, there are 3,500 hit pixels in a 
neutrino event. 

2. Finding Pixel Clusters A pixel cluster is determined by a group of 
neighboring hit pixels. At this stage, isolated single pixel hits, coming from 
random electrical noise of CCD, are rejected. 
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Fig. 4. An Npix distribution with the corrections apphed for both attenuation and 
dependence on the incident angle as given in Eq.(l). 

3. Identifying Hit Fibers A hit fiber is identified by requiring at least 
one or two hit pixels of cluster within the fiber. The hit fiber position in the 
CCD coordinate is determined from a mapping table that is obtained from 
the EL calibration. 

4. Making Hit Fiber Clusters A fiber cluster is determined by a group 
of neighboring hit fibers. 

We apply two kinds of filtering algorithms to the fiber clusters to improve the 
signal-to-noise ratio: 

• We require that a fiber cluster should consist of hit fibers in both upper 
and lower layers of a fiber sheet. This eliminates fiber clusters of small pixel 
images which come from dark current in an IIT. 

• When a big pixel cluster lies on the adjacent fiber sheets in a bundle and 
makes multiple fiber clusters, we choose only fiber clusters containing peak 
of ADC brightness. This removes fake fiber-clusters introduced by the dense 
arrangement of fiber sheets in a bundle. 

We consider a fiber cluster passing through these two filters as a SciFi "hit" . 
Every hit position in the CCD coordinate is then converted into the position 
in the detector coordinates or y,z). The hit position is defined as the 
centroid of a fiber cluster, weighted by the ADC value of the hit pixel. 
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3.2 Number of pixels per hit(Npix) and its corrections 



A found hit gives the number of pixels for the parent pixel cluster, Npi^. LED 
calibration data show that Npi^ has good linearity with respect to the number 
of photons injected to IIT. Since the intensity of observed light is proportional 
to the energy loss in the fiber, Npt^ is expected to be a good estimator of the 
energy of the parent particle. 

Fig. 4 shows an Npi^ distribution for cosmic-ray muons. Based on a recon- 
structed track, Npix can be corrected to account for both the attenuation of 
light in the fiber, and the dependence of light yields on the incident angle (6'j„) 
of a particle on a fiber sheet: 



-^pixcorr 17/ J\ ^ C0s9in X Npix- (1) 

The first term is the correction for attenuation of the scintillation light in 
the fiber and reflection at the far end of fiber coated with aluminum [5]. The 
correction is normalized by the value in case of the light from the detector 
center (/o=130cm). Thus, F{1) is given by: 



F(l) = expi--^) + R ■ ea:p(-^-^^^), (2) 

where / is the distance from the hit position to the near end of a fiber, 
lg{=100cm) is light guide length from the near end to the IIT, and A(=323cm) 
is the attenuation length. i?(=0.74) is the reflection coefficient at the far end 
of a fiber, and L(=254cm) is the length of a fiber. 

The second term in Eq. (1) is the correction for the dependence of light yields 
on the incident angle with respect to the fiber sheet. Here, 6in=0 is perpen- 
dicular to the fiber sheet. Npi^corr is thus equivalent to Npi^ when a particle 
passes fiber sheet at the center of the detector x/y planes with 6*^=0. 

Figs. 5(a) and (b) show average values of Npi^ {{Npi^)) as a function of hit 
position along the fiber and as a function of l/cos{9in), respectively. 

^ We employ the following Poisson function to fit the Npix distribution: 

F{Npix; A, a, {Npix)) = A ■ ^ , /i = {Npix) /a, and x = Npix/a, 

where A is an overall normalization factor and a is a conversion factor. (Npix) is 
the average value of the Npix distribution. 
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Fig. 5. An average value of Npi^, {Npi^), as a function of (a) hit position along the 
fiber and (b) l/cos6in, where 6in is the incident angle with respect to the fiber sheet, 
respectively. Solid (open) circles are before (after) the correction. The negative x, y 
coordinates in (a) correspond to positions close to readout end of the fiber, and a 
larger value of {Npi^) is expected before the correction. 

3. 3 Simulation of CCD hit pixels 



To reproduce SciFi hits accurately in the K2K Monte Carlo simulation, we 
include simulation of hit pixels on CCD due to charged particles and noise in 
the following three steps: 

1. Charged particles are tracked through the detector using GEANT[6]. The 
energy loss of a particle traversing a fiber is computed and then converted 
to the number of photoelectrons observed by the IIT/CCD imaging system, 
taking into account attentuation along the fiber, reflection by the aluminium 
coated at the end of a fiber, and the quantum efficiency of the IIT photo- 
cathodes. 

2. The photoelectrons are converted to an image on the CCD camera surface. 
Individual photoelectrons are distributed around the fiber center according 
to a gaussian distribution with parameters determined using data of a sin- 
gle photoelectrons from LED calibration run. The size of the pixel cluster 
corresponding to the single photoelectron is then chosen from distributions 
of the LED data, and the pixel cluster is simulated accordingly. The energy 
scale is tuned for each IIT using cosmic-ray muon data. 

3. Random noise in a pixel cluster is generated for each IIT based on neutrino- 
beam data. 

Two kinds of noise pixels are generated. One is isolated single hit pixels around 
the main pixel cluster of a charged track in step 2.. The other is random noise 
pixels in a pixel cluster generated in step 3.. The total number of noise hits 
averaged over IITs is 1.67 for MC and 1.72 for cosmic-ray data. Based on the 
pixel simulation, the average number of noise hits is expected to be 70 per 
event. In the cosmic-ray data we find 52 hits per event due to the random 
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Fig. 6. Hit finding efficiency Chu obtained from cosmic-ray data 

noise of CCD, and 17 hits per event due to the beam associated background. 
The noise rate is quite low, considering the total number of fibers (274,080). 



3.4 Efficiency and Stability of Hit Reconstruction 



We can obtain the hit finding efficiency ehu which is defined as: 



Number of observed hits 
Number of expected hits' 

where the expected number of hits (denominator) is obtained from the infor- 
mation of reconstructed tracks. Fig. 6 shows ehu for each IIT. Data and Monte 
Carlo shows a good agreement. 

Cosmic-ray data are taken during the normal data-taking period, roughly 
once every two weeks. The efficiency is monitored using the data as a check 
of detector stability. 

We have performed a similar check on the stability of hit finding efficiency 
using the neutrino data. Clean single track events are selected by requiring 
a track with the number of hit layers larger than nine, and its matching to 
activity in the downstream VETO counter. 

Fig. 7 shows hit finding efficiency, which is obtained by taking average over all 
IIT's, as a function of data-taking period. The dotted line at 96.8% represents 
an MC prediction. For the cosmic-ray data, we find an average value of e^it to 
be 96.5 ± 0.5% for all IIT's and all time periods. 0. 

^ Tliis is in good agreement with a simple estimation. When a minimum ionizing 



9 



o Beam 

• Cosmic-ray 



ON © © 

Fig. 7. Hit finding efficiency as a function of data taking period for the neutrino 
data (open circle) and cosmic ray data (solid circel). The dotted line represents an 
MC predicted hit finding efficiency. 

4 Reconstruction of Track and Vertex 

4-1 Track Reconstruction Algorithm 

The track finding algorithm in the SciFi is optimized for the charged cur- 
rent interactions which contain a muon in the final state. The track finding 
algorithm consists of three steps. 

First, the track finder searches for 2-dimensional (2D) track candidates in 
the {x, z) or (y, z) planes which traverse at least three layers of SciFi. A 3- 
dimensional track is reconstructed by selecting the best combination of two 
2D tracks. In order to reduce the combinatorial fake tracks, the starting and 
ending hit layers are required not to be different by more than one layer, and 
the overlapping of the two 2D tracks is required to be at least two layers. 

Secondly, the combinatorial background or noise tracks are removed using 
other components of the near detector (VETO and MRD). A track candidate 
in the SciFi is required to match with both a hit in the VETO and a track in 



particle crossing the fiber sheet creates ~ Sp.e. {4p.e. each for upper and lower lay- 
ers), hit coincidence probability between upper and lower layers is (1. — e~^)^=0.96. 
The requirement of the coincidence between two fiber layers is chosen to keep the 
noise hits at the level of 70 hits per event. 
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the MRD. 



Once a good primary track is found, the most upstream hit-position of the 
primary track is considered as a neutrino interaction point. The track finder 
starts searching for extra tracks near the neutrino interaction point. These 
tracks are mostly short tracks with two or more SciFi hits. This method of 
extra track finding turns out to be an effective way to find short tracks of 
protons and/or pions from the neutrino interactions. AU the reconstructed 
tracks are subjected to a "global fit" , a straight hne fit to the found hits. 

4.2 Alignment and Position Resolution 

The geometry of the SciFi and its support structure was surveyed by mechan- 
ical means after their installation was completed. Based on the mechanical 
survey, the alignment of fibers is estimated to be good to within 5mm. 

A better alignment was obtained by minimizing the sum of the squared hit- 
residuals using the global fit for cosmic ray tracks. From the mean shift of the 
hit residuals, we find the accuracy of the alignment to be ~400/xm. 

Multiple scatterings occurring in the SciFi is not negligible compared with the 
size of a fiber. A Kalman filtering technique [7] [8] [9] can be employed to take 
the multiple scattering into account and, in principle, removes the multiple 
scattering effect during track fitting. The technique calculates the error of a 
track fit relatively quickly, and avoids the tedious matrix inversion by the 
global fit. A track fit using the Kalman filtering technique is used to estimate 
the detector resolution more accurately and to obtain a better alignment. The 
implementation adopts the prescriptions described in Refs. [9] and [10]. 

Some curvature was introduced in the fiber sheets due to bending when they 
were glued to a fiber module. Production of the fiber modules is described in 
Ref. [5]. The curvature of a fiber sheet causes displacement from the aligned 
fiber-position and results in a rotational misalignment. The bending of every 
fiber sheet was measured after each fiber module was assembled. The curvature 
of a fiber sheet was then obtained by fitting the measured fiber bendings to a 
forth-order polynomial, and used to correct the rotational misalignment. 

After the Kalman filtering and the curvature correction, the accuracy of align- 
ment becomes better than 50/xm, estimated from the mean residual shift. Fig. 8 
shows distributions of the mean and standard deviation of the hit residuals 
for the mechanical survey alone, the alignment using the global fit, and the 
alignment using the Kalman filtering fit. 

The mean of hit residual after the mechanical survey alone shows a variation 
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Fig. 8. (a) Means of hit residual and (b) standard deviations as a function of layers 
in y — z plane. Points are for the mechanical survey alone (open star), the alignment 
using the global fit (open circle), and the alignment using the Kalman filtering (solid 
circle) . 

as large as 1cm. After the alignment using the global fit, the variation becomes 
smaller than 400//m. After the alignment using the Kalman filtering fit, the 

variation becomes at most 50;um. The standard deviations of the hit residual 
also exhibit such improvements of fiber alignment as the mean values. There 
are a few layers with large standard deviations of hit residual, compared with 
other layers. These are the outer most layers and indeed bent more than the 
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rest of inner layers. 

The alignment is checked using periodically taken cosmic-ray data. Fig. 9 
shows the variation in the mean values of hit residuals since May 1999, and 
the variation is found to be less than 100//m. 
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Fig. 9. The variation in the mean of hit residual, obtained from cosmic-ray data, as 
a function of layers in x — z plane (1-20) and y ~ z plane (21-40). The variation is 
taken relative to May 1999. The dotted line represents the variations of itlOO/um. 



The neutrino data with the Kalman filtering fit are used to obtain the best 
estimation for the position resolution. Fig. 10 shows the hit residual distri- 
bution using neutrino data. The position resolution is 0.73mm (0.61mm) for 
muon energy < IGcV ( > IGeV). The overall position resolution is 
0.64±0.07mm. Fig. 11 shows the energy dependence of the resolution. The 
effect of multiple scattering on the position resolution decreases as the muon 
energy increases. The MC prediction explains the overall energy dependence. 
The difference between data and MC, ~100/xmm, is regarded as the system- 
atic error of our detector. The position resolution of 0.6mm is regarded as that 
of the SciFi without the multiple scattering effect 
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Fig. 10. The hit residual distributions using the neutrino beam data with the Kalman 
filtering fit. The curve is a single Gaussian fit to the data with cr = 0.73mm (0.61mm) 
for muon energy < IGeV ( > IGeV). 

4-3 Track Finding Efficiency 

The efficiency for finding a single muon track is estimated from cosmic-ray 
and MC events. The cosmic-ray events are selected by requiring associated 
activities in both upstream and downstream VETO counters. Furthermore, 
tracks are required to have at least 15 SciFi hit layers from the upstream to 
the downstream. In order to study dependence of the track finding efficiency 
on track length, we artificially eliminate hits in the upstream side of the track. 
For example, if it is desired to estimate the efficiency of the track with ten 
SciFi hit layers using 19-layer long track, the first to the nineth hit layers are 
eliminated before applying a track finding algorithm. For MC events, a muon 
track is required to have at least three hit-layers within the SciFi fiducial 
volume, and tested by a track finding algorithm. 

Fig. 12 shows a comparison of single track finding efficiency as a function of the 
track length in unit of number of SciFi layers for the cosmic ray data. Track 
finding inefficiencies come from failures in the hit ffnding and imperfection of 
the track reconstruction. The track finding efficiency is ~70% (~87%) for the 
track length of three (four) hit SciFi layers and becomes close to 100% for 



14 



1 

0.9 

0.8 



03 



3.6 



:g 0.5 

W 
0) 

DC 0.4 

T. 0.3 

0.2 
0.1 





I v„ beam data 
MC 







E, (GeV) 



Fig. 11. The energy dependence of position resolution. 

longer tracks. 

An average track finding efficiency for the cosmic ray events is estimated by 
weighting the track length distribution obtained from the MC. We find the 
average track finding efficiency to be 93.6±0.8% for cosmic-ray tracks with 
three hits or more. The average track finding efficiency is estimated to be 
92.7±0.5% for the MC events. 



4-4 Vertex Reconstruction 



An accurate fiducial volume is crucial for measuring the cross section of neu- 
trino in the SciFi. Therefore, the event vertex has to be determined precisely. 

The vertex position is determined by the reconstructed tracks. The event 
vertex for a single track event is found by simply taking the mid-point of 
the upstream water container nearest to the first hit in the z-direction of the 
reconstructed track. The event vertex for a multiple track event is determined 
by taking average of the z positions at intersections of a primary track and 
the rest of shorter tracks. Fig. 13 shows a comparison of x- and ;2-vertex 
distributions between data and MC for both single and multiple track events. 
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Fig. 12. The primary track finding efficiency as a function of track length using 

cosmic ray data. SoHd rectangle, open circle, and solid circle are for 1999, 2000, and 
2001 cosmic ray data, respectively. The track length is given in units of number of 
SciFi layers. 



Events with tracks ranged out at the MRD are used for the comparison and 
the comparison shows a good agreement between data and MC. 

For two track events, vertex difference between {x, z) and (y, z) planes can 
provide a measure of goodness for the vertex reconstruction. Fig. 14 shows 
a distribution of 2;-vertex difference between vertices found in the (x, z) and 
(y, z) planes for events with multiple tracks. The event z-vertex is found by 
taking the average of the z positions determined for all such combinations of 
tracks. The event vertex in the x and y coordinates is calculated based on the 
event ^-vertex and the primary track. 

The vertex resolution is estimated by two track MC events of neutrino charged- 
current interactions. The vertex resolution is degraded mainly due to track 
mis-reconstruction. Fig. 15 shows the difference between the measured vertex 
position and the true vertex position in y, and z for the two track events. 
The estimated vertex resolution for the two track events is found to be 0.6 cm 
in z and 0.4 cm in x or y. 
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Fig. 13. A comparison of x- and z-vertex distributions between data and MC for 
both single and multiple track events. 



4-5 Track Counting 



An accurate reconstruction of the track multipUcity of an event is crucial 
for the measurement of neutrino energy spectrum. Charged current quasi- 
elastic events have especially simple kinematics which can be used to determine 
the incoming neutrino energy. These events always have one or two tracks 
depending on whether the outgoing proton is detected. 

The efficiency of finding multiple tracks is studied using the MC CCqe events, 
and defined to be the fraction of reconstructed events relative to the generated 
events for a given track multiplicity. Both the reconstructed track and the 
generated track are required to have three or more hits. The efficiency is 
estimated to be 93±2% and 86±4% for single track CCqe events and two 
track CCqe events, respectively. For single track events, the uncertainty in the 
track counting efficiency is dominated by the systematic error due to the gain 
variation of the IIT's, and estimated to be 1%. 

For two track events, the sources of the uncertainty in the track counting ef- 
ficiency are difference in noise between MC and data, tracks due to noise hits 
around the vertex, and noise hits near the primary track. Due to the differ- 
ence in noise between data and MC, the track reconstruction may perform 
differently. 
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Fig. 14. The z-vertex difference between vertices found in the (x, z) and {y, z) planes 
for two track events. At least three hits are required for both primary track and 
second track. Points are for data and the histogram is for MC. 

In order to estimate the systematic uncertainty due to difference in noise 
contamination between the data and MC, noise contamination is varied in 
the MC by ±50% from the default value. After the variation, track finding is 
performed again. We find the track counting efficiency varies by at most 1.8% 
with respect to the default noise containmination in the MC for both single 
track and two track events. 

Single track events can sometimes be reconstructed as two track events be- 
cause of noise hits around the vertex. The size of this effect is studied using 
MC generated single track events with three or more hits. The probability is 
calculated by counting number of reconstructed two track events, and found 
to be 2.7±0.1% for the MC CCqe events. 

On the other hand, if there are nearby noise hits along the primary track, 
single track events might be reconstructed as two track events. This type of 
nearby noise hits contributes to the uncertainty in track counting for two track 
events. For single track events, the vertex position is relocated to the middle 
of the track and the track finding is repeated around the new vertex. The 
uncertainty in track finding due to accidental noise hits around the primary 
track is found to be less than 1%. 
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Fig. 15. The vertex resolution of two track charged-current events. 

The track counting inefficiency can be explained by loss of tracks or additional 
fake tracks through reconstruction. The efficiency of finding the second track 
is calculated as a function of number of hits in the SciFi layers using MC. 
Events with two generated tracks are reconstructed, and only those with a 
primary track and a vertex are selected. The second generated track should 
have at least two hits. The efficiency of finding the second track is defined to 
be the ratio of number of events with two reconstructed tracks to the number 
of events with two generated tracks. 

Fig. 16 shows the efficiency of finding the second track as a function of number 
of hits in the SciFi layers. The upper plot is for all types of MC events, and 
the lower plot is for the CCqe events. 

For the second tracks with three or more hits, the inefficiency comes from 
hit finding inefficiency and limited performance of the track reconstruction as 
for the primary track. For the second tracks with two hits, the inefficiency 
is mainly due to noise hits around the vertex or near-by noise hits along the 
primary track. An average efficiency of finding the second track with three or 
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Fig. 16. The efficiency of finding the second track as a function of number of hits 
in the SciFi layers. The upper plot is for the all events and the lower plot is for the 
CCqe events. 

more hits is found to be 82.1±0.4% for all types of MC events and 90.1±0.5% 
for MC CCqe events. 



5 Event Reconstruction 



In the previous section, we described the reconstruction of tracks and ver- 
tices. In this section, we demonstrate their performances using distributions 
of typical physics variables. 

Fig. 17 shows a typical i/^ CCqe candidate event. The long track is the muon 
which leaves hits in the VETO and LG, and ranges out at the MRD. The 
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Table 1 



Summary of tracking performance of the SciFi. Position resolution is measured using 
the neutrino data. Vertex resolution is estimated using MC events. Track counting 
efficiency is calculated for the CCqe events. Errors on the hit efficiency, tracking 
efficiency, and track counting efficiency are the quadratic sums of statistical and 
systematic errors. 

short track is most likely a proton. 

Most of the protons from neutrino interactions have relatively low momenta 
and thus make short second-tracks in the SciFi. Fig. 18 shows the track length 
distribution for the second tracks in the 2-track sample. Both Monte Carlo 
prediction and data agree well, demonstrating the validity of the neutrino 
event generation and the tracking method. 



6 Conclusions 

The K2K long-baseline neutrino experiment has taken data for three years. 
The SciFi is used to identify charged particle tracks in the K2K near detector. 

The SciFi hits are calibrated using cosmic ray data and electro-luminescent 
calibration data. MC simulation reproduces the distribution of the number of 
hit pixels in the cosmic ray data very well. The stability of the SciFi is checked 
periodically by obtaining the distribution of hit efficiency and calculating the 
mean of the hit pixel distribution. The stability is better than 5% and 1.2% 
in the mean of pixel distribution and the hit efficiency, respectively. 
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Fig. 17. A candidate i^^ CCqe event. The long track is the muon and the short track 
is most hkely a proton. 



A track reconstruction algorithm, optimized for charged current neutrino in- 
teractions, is employed to find tracks and event vertices. The track finding 
efficiency is monitored and found to be stable throughout the data-taking 
period. Using the estimated vertex resolution, the vertex reconstruction alo- 
girthm is found to be good enough for defining the event fiducial volume with 
an accuracy better than 1%. A track fit using a Kalman filtering technique is 
used for the best estimation of hit resolution, and found to improve the fiber 
alignment. 

Accurate track counting is crucial for the measurement of neutrino energy 
spectrum and in achieving the physics goals of the K2K experiment. The 
track counting uncertainty is found to be at 2% and 4% level for the single 
track and two track events, respectively. Overall tracking performance of the 
SciFi is summarized in Table. 1. 

The performance of the SciFi tracker is well demonstrated by physics variables 
of neutrino interaction. The distributions of variables for muon track and 
proton track show a good agreement between data and MC events. 



22 




20 40 60 80 100 120 140 160 ISO 200 



Track Length (second track) (cm) 

Fig. 18. A comparion of track length for the proton between neutrino data and MC 
events in two track sample. 
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